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Unique continuation property of solutions to 
general second order elliptic systems 

N. Honda* C.-L. LiU G. Nakamura^ S. Sasayama^ 
Dedicated to the memory of Prof. Kenjiro Okubo 


Abstract 

This paper concerns about the weak unique continuation property of solu¬ 
tions of a general system of differential equation/inequality with a second order 
strongly elliptic system as its leading part. We put not only some natural as¬ 
sumption which we call basic assumptions, but also some technical assumptions 
which we call further assumptions. It is shown as usual by first applying the 
Holmgren transform to this inequality and then establishing a Carleman esti¬ 
mate for the leading part of the transformed inequality. The Carleman estimate 
given via a partition of unity and Carleman estimate for the operator with con¬ 
stant coefficients obtained by freezing the coefficients of the transformed leading 
part at a point. A little more details about this are as follows. Factorize this 
operator with constant coefficients into two first order differential operators. 
Conjugate each factor by a Carleman weight and derive an estimate which is 
uniform with respect to the point at which we froze the coefficients for each 
conjugated factor by constructing a parametrix for its adjoint operator. 


1 Introduction 

Let A G N and be a domain in with n > 2. Snppose a C" valned vector 
fnnction u{x) = ■ ■ ■ , un { x )) G G^(r2) satisfies the differential ineqnality 

\Lu{x)\ < c ^ \d''u{x)\ in hi for some c > 0, (1.1) 
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where L denotes the N x N system 


n N 

{Lu)a{x) = EE Ci^p{x)djdiuii{x) (1 < a < iV) with dj = dx^ (1 < j < n) 

j,l=l / 3=1 

and d = (9i, • • • , 9„). Through out this paper, any vector is considered as a column 
vector. We first assume the following basic assumptions which are the smoothness, 
symmetry and strong ellipticity conditions on L. 

Basic assumptions on L 

(i) smoothness condition : Each is a real valued Lipschitz continuous function 
in n. 

(ii) symmetry condition : 

= Cpiix) {x e n, 1 < j, i < n, 1 < a, 13 < N). (1.2) 

(iii) strong ellipticity condition ; There exists 5 > 0 such that for any vectors a = 
(oi, • • • , Oat) G and 6 = (6i, • • • , 6„) G M"' 

C^Jp{x)aabjaybi > 6\a\‘^\b\‘^ {x G Q). (1.3) 

1 < q ;, i<n 


The aim of this paper is to consider the weak unique continuation property simply 
abbreviated by UCP of fll.ll) . i.e., if u{x) is any solution of fll.ip which vanishes in a 
non-void open subset of hi, then u{x) vanishes identically in hi. We prove the UCP for 
the general differential inequality fll.ll) with L satisfying the above basic assumptions 
and also some technical assumptions which we call further assumptions given in the 
next section. Unlike the case for scalar partial differential equations/inequalities, 
there are very few results known about the UCP for systems of partial differential 
equations/inequalities. 

A very important example of a strongly elliptic system of partial differential equa¬ 
tions satisfying the above three conditions is anisotropic elastic system of partial 
differential equations which arises in mechanics of materials and geophysics. In this 
case n = N and is the elastic tensor held. When the medium is isotropic, the 
UCP has been established in n, i. On the other hand, for the anisotropic medium, 
there is a rather general result known for the two dimensional case ( 0 ). but there 
are only very few results known for the three dimensional case. 

A general and powerful method of proving UCP was pioneered by Calderon ([2]) 
and has been generalized by several authors ([T3| and references therein). But this 
is basically for scalar partial differential equations/inequalities and in order to apply 
this method to systems of partial differential equations/inequalities one needs to 
diagonalize the leading part of system of partial differential operators to a system 
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of pseudo-differential operators whose characteristic roots of the principal symbol 
are smooth and they are simple for real characteristic roots and at most double for 
complex characteristic roots. But these conditions are usually not satished for general 
anisotropic system of elastic partial differential equations and not even for general 
transversally isotropic elastic systems. 

Our method, as in previous results for UCP, relies on a suitable Carleman esti¬ 
mate. We adapted the idea of Vessella et ah ([12]) to reduce the derivation of a 
Carleman estimate to that for the operator with constant coefficients obtained by 
freezing the coefficients of the leading part of fll.ip undergone the Holmgren trans¬ 
form. Then, we adapted the idea of Sogge m) to derive a Carleman estimate for 
this operator with constant coefficients which is uniform with respect to the point at 
which we froze the coefficients. More precisely, we hrst factorize this operator with 
constant coefficients into two hrst order differential operators. Next we conjugate 
each factor by a Carleman weight with large parameter k. A similar factorization 
was used in [8] for a general elastic system of partial differential equation to study a 
layer stripping method for its associated inverse boundary value problem. Finally, for 
each conjugated factor we construct a suitable parameterix to derive the Carleman 
estimate. 

The rest of this paper is organized as follows. Section 2 gives some preliminaries 
necessary to state our main results. Especially the aforementioned further assump¬ 
tions is given there. In Section 3, we state a Carleman estimate and prove UCP using 
this estimate. The Carleman estimate will then be derived in the remaining sections 
based on the argument outlined in the previous paragraph. More precisely, as fol¬ 
lows. Section 4 provides the factorization of constant coefficients system of partial 
differential operators L obtained by freezing the coefficients of L which undergone the 
Holmgren transform into the product of two factors which are the hrst order differen¬ 
tial operators in t and hrst order pseudo-diherential operators in the other variables. 
Also, they do not depend on t and other variables. Here t is the variable whose axis 
gives the direction to which we want to do UCP. Then these factors are conjugated 
with a Carleman weight Wkit) and obtain two pseudo-diherential operators Pk,g and 
Pk^b with large parameters k which do not depend on t and other variables. Pk,g 
is good because it is elliptic, but Pk,b is bad because it is not elliptic. The right 
parametrix Sk of (T — is given in Section 5. The error of Sk is estimated in 

Section 6 to show that the fundamental solution of (T — t)~^P^^ can have the same 
estimate as that of Sk- The estimate of Sk will be gi ven in Sections 7, 8. Based on 
these estimate for Sk and that of we give a Carleman estimate for L in Section 
9. Finally in Section 10, we give a Carleman estimate for L using the partition of 
unity. 
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2 Preliminaries 


In order to state our main result in the next section we will give some further assump¬ 
tions on the operator L. Let L be a hypersurface in R” intersecting with Take 
Xq € r and its open neighborhood U G fl. Denote the symbol of —L by M = M{x,()- 
For a unit conormal vector 77 of F and any ^ G R"" independent to 77 , let C = ^ 

with A G R and write M in the form 

M = TX‘^ + AX + Q, (2.1) 


where 

n 

T = T{x,ri) = C^^p{x)r]j 77^; a I 1, • • • , ^ 1 • • • , iV), 

i/=i 

A = A(x, 77,0 = h, 0 + h, 0, 

n 

R = R{x, r],0={Y a ; 1, • • • , iV, ^ 1 ■ ■ • , iV) , 

n 

Q = Q{x,o = {Y « ; 1 , • • •, iv, ^ 1 • • •, AT) , 

i,^=i 

where the superscript ”T” denotes the transpose of matrices. 

It is easy to see that the positivity of T follows from the strong ellipticity condition 
of fll.31) . Hence and exist and we can consider 

H{X) = H{x, X, 77,0 = T-V2MT-V2 = + H,{x, 77, OA + H,{x, 77,0 (2.2) 

with 

Hi = t-^RaT-^I‘^ and H 2 = (2.3) 

It is known that H{X) has a unique factorization 

H{X) = {X-S\){X-Si) (2.4) 

such that the spectrum Spec Si of Si satishes Spec Si C C+ := {z G C : Im 2; > 0}, 
and this unique Si can be given by 

Si = Si{x,v,o--=i{ (Hicruoif H{c)-^do-\ 

J7+ J7+ 

where 7+ C C+ :={CgC: Im(C>0}isa closed contour enclosing all the roots of 
det(iL(a:, A, 77,^)) = 0 in A and S']" denotes the adjoint of Si ([!]). Further, 

Z = Z{x, 7], 0 = Im(S'i) (2.5) 
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is positive-definite ([5]) and for V := Re(S'i), it can be easily proved from S^Si = H 2 ., 

^l/2y^-l/2 

is symmetric. Hence by 

( 2 . 6 ) 

'H = 'H(A) has a mnch more nice factorization: 

n^X) = {\-B*)Z-\\-B) (2.7) 

with the complex symmetric matrix 

B = B{x, 7],C! = BR + iBi, (2.8) 


where Br^ Bj are given by 

Br = Br{x,v,0 = Re(Zi/2^iZ-i/2), Bj = Bi{x,r],0 = Z. (2.9) 
Clearly the spectrums Spec(i?/j ± iBj) of Br ± iBj satisfy 

Spec(i?ij ± iBj) C C± := { 2 ; G C : ± Im ^ > 0}. (2-10) 

Now let Ph = Ph{x,7],^) be the projection of B given by 

F, = P,{x, V, 0 = (2vrz)-' f {(- B)-^dC (2.11) 

•Jih 

with a positively oriented small circle 7/1 = 7 h(x, 7 ,^) in C+ enclosing the eigenvalne 
Xh = Xh{x, h) 0 ^ of B but excluding other eigenvalues Xk = Xk{x, h) 0 ^ 7 ^ 

h, 1 <k < s). Here the multiplicity ruh of Xh satishes rrti = N . Then the 

further assumptions we put on L are as follows. 

Further assumptions on L 

(i) Each Ph{x, rj, is uniformly bounded for a.e. x E U and all ^ G M”\{0} linearly 
independent to 7 . 

(ii) B{x, rj, is diagonalizable for a.e. x E U and all G \ {0} linearly indepen¬ 
dent to rj. 


Remark 2.1 

(a) By Theorem 1.9 in m, for each almost everywhere fixed x E U, the assumption 
(i) implies that all the eigenvalues Xk{x,ri,f) {1 < k < s) of B{x,ri,f) have to be 
analytic in f G M" \ {0} which is linearly independent to rj. Needless to say this is 
equivalent to saying the same for all the eigenvalues of M. 
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(b) The assumption (ii) can be given in terms of the principal symbol M of —L. That 
is for a.e.x E U, every G M” \ {0} linearly independent to rj satisfies 


dime Ker M(a;, At] + .^) = 2iV, (2-12) 

AeA(a:,^) 


where 

A{x,f) = {A G C : det M{x, Xp + ^) = 0} (2-13) 

and dime means the dimension for complex vector spaces. 

(c) Suppose we have the factorization fl2.4p for fl2.2|) . Then note that the symmetry 
of Y is equivalent to the commutativity of H 1 H 2 = H 2 H 1 . We refer this as com¬ 
mutativity. If this commutativity holds for any x E U and ^ E \ {0} linearly 
independent to fixed p, there exist an orthogonal matrix G = G{x, p) such that 


H 2 = 


where Si and S 2 are diagonal matrices. Then by the strong ellipticity condition, we 
have Til ~ < 0- Thus, we can factorize H itself very nicely as follows. 


H = 


■yT 


A" + AGSiG' +GT2G' 

G(A2 + ASi + S2)G^ 

G(A + f + El + 


(A - B*){X - B) with B = -G(f 




It is easy to check that the complex symmetric matrix B = Bji + iBj satisfies the 
further assumptions. Furthermore, if we can assume G{x, f, p) is smooth in x E U 
and G M"' \ {0} linearly independent to p, we can transform the operator L to a 
system of pseudo-differential operators of order 1 which is a differential operator in 
xi and has diagonal principal part. By applying Calderon’s argumentwe can 
have a Carleman estimate which implies the uniqueness of the Cauchy problem for the 
initial hyperplane F. If we have such a situation for every hyperplane T intersecting 
we have UCP in fl. However, this smoothness condition on G does not always 
hold. Example 4 given below shows such a case. 


Next we will give several examples which satisfy all the basic assumptions and 
further assumptions. For simplicity we take p = (1, 0, • • • ,0) and write A = ^ 1 , ^ = 

(6, • • • ,^n)- 

Example 1. Let 


V 


The eigenvalnes of B are 




7^2-V3\f2\ , 3|e2| + \/36. re2 + A/3|e2| , 3161-\/3e 


+ 


+ 
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and the corresponding eigenvectors are 


(10161 + 26h 5^36 - 5|6I + (5^161 - 116)*), 
(-10161 - 26*, 5\/36 + 5|6I + (563161 + 116)*)- 


The associated H is 


H 



-lo) + 






2 

2 - 


and it does not have the commutativity. 

Example 2. Perturb the previous example as follows. Let 



with El, e 2 i q & ^ satisfying \q H-| > 5 > 0 and \q -> 5 > 0 for some 

The eigenvalues of B are 


A± 

146 Tc|6l * 

4 ^2(1 +a) 


1^(3 + iqei + q^el - 2qe2 - q^el)\^2\ ± 


^(\/& + 6 + 2 a) 


where 


a — qi^^i — ^ 2 ), 5 — a T 4a^ T 40a T 72a T 36, c 



— — 2a. 


Corresponding eigenvectors are 


(/3, Aj- o), 


where 


_ {(P + d{2qei + 1 ) + qsi — 1)6 + *|6 

d(3 + q(y£\ + £ 2 ) + 2p) 

{AP + q{ei + £ 2))6 + *(2 + 5'£i)|6| 

(y = - 

d{3 + q{ei + £ 2 ) + 2d?) 
d = 61 T q^i T ‘?q^2 + q‘^^1^2- 


The associated H is 


H 


/—4 0 f flQ + Aq£i + q^£\ Q + q£i + q£2 \ ^2 

yO -lOy'^^'^^ y Q + q£i + q £2 28 + 2ge2 + 6^2/ 


This H also does not have the commutativity. 
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Example 3. Let 


/I 0 0\ /I 0 0 \ 

F = 0 3 16+0-1 16, 

\0 -1 0/ \0 -1 -4/ 



The eigenvalues of B are 


6 + 'Cs + i\l3^2 + 26^3 + Cl, 

^ 36 - 56 ± & ^ , 9(36 + 56^1 + 3^16 + W ± ± 566 

2 6(6 + ■ Cs ) \/36 + 2 ^ 2^3 + .^1 

where 


b := a := 526^ + 726^6 + 726^3 + 366^1 + 96"- 

The corresponding eigenvectors are 


(1,0,0), (/9, A±-a), 


where 

13 '■= + 26^3 + ^3, “ := —^^ + ‘^i\J3>^2 + 26'C3 + ■C; 

The associated H is 

/-26 - 26 0 0 \ 

H=ej + ^i\ 0 -66 + 26 0 

V 0 0 86/ 

/4e| + 46e3 + 6a 0 0 \ 

+ 0 22^1 + 466 + eel eel + 86e3 + 4e| . 

V 0 6ei + 8e2e3 + 4ei 4ei + 4e2e3 + i8e37 

This H does not have the commutativity. 

Example 4. Let 

// = (6 - - B) with B = Br + Wj, 

where 
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(2.14) 



The eigenvalues of B are 


|-\/(2^26)^ + (a^sCs)^ + i\lil + Ci, 'J + (a^sCs)^ + + C|, + ^|| , 

and the corresponding eigenvectors outside the set {x2 = 0:3 = 0} are given by 

Vl := (x2^2, -\/(x2^2)^ + (o^sCs)^) , 

V2 ■= (3^26, 2:3^3, a/(0:26)^ + (2^36)^) 

and 

t'a := (3:36, -X2^2, 0), 

respectively. This example shows that all the conditions are satished if 0:2X3 7^ 0 and it 
does have the commutativity. However, there is no continuous diagonalization where 
we have X2 = X3 = 0. As a matter of fact, the 1-dimensional vector bundle spanned 
by the vector V3 on {x‘l + x\ 7^ 0} never extends to the whole space continuously since 
the point (x3^3, —X2C2) turns once around the origin in if ,^2^3 7^ 0 and (x2, X3) 
moves once around the origin. 

Example 5 . Let 

H= [iiJ-(A-^I^EAy^y+ iAs/^TWA 
X [?a - A-'/‘^EA'l^ - iA^(f+(l] 

= - (A-V2^H1/2)T _ 

+{A-^/^EA^/^y{A-^/^EA^/^) + (^1 + ^|)A2, 

/I 0 \ /I 0 0 \ 

where A = 0 1 -|- (x2X3)^ 0 and E = 0 X2^2 . We assume 

yO 0 1 - {x2X3yJ \0 X3^3 -X26/ 

that, in what follows, (x2, X3) is sufficiently close to the origin. The eigenvalues of 
Si = Y + iZ = which is in fact equal io B = Bj + iBj are 

given by ir^, 

ir^ — r|x|x| + 2i^2r^X2xl + 

and 

ir^ -f ^J-r'^xjxj + 2 i^ 2 r^xlxl + 

where we set := and := ^/(x2^2y~+Jx^^^■ The corresponding 

eigenvectors on X2X3^2^3 7^ 0 are Vi := (1, 0, 0), 

^2 := ( 0 , -sj-rlx^xl + 21^2 r^x^xl + rl^ + i r^xjxl + ^ 2 X 2 , ^ 3 X 3 ) 

and 

^3 := (0, sj-rlxlx'^^ + 2 i ^2 r^xlxl + rl^ + i r^xlxl + ^2X2, C3X3), 
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repsectively. This example shows that all the conditions are satished if X2X3 7^ 0 
and it does not have the commutativity. Also, there is no continuous diagonalization 
where we have X2 = X3 = 0 . These observations follow from the fact that there exists 
a constant M > 0 such that, on ^2 + ^3 = 1 ) have 


1^2 - (0, + ^23^2, ^32^3)1 < Mrl ^ 


and 


1^3 - ( 0 , + ^2X2, 63^3)1 < Mrl^. 


Through out this paper we will always assume the above basic assumptions and 
further assumptions. 

3 Main results 

First of all as it has been already mentioned before in the last remark given in the 
previous section, the UCP in ^2 follows from the following uniqueness of the Cauchy 
problem (see [H]) for any hyperplane intersecting with ^2. Also, the uniqueness of 
the Cauchy problem follows from a suitable Carleman estimate given below in this 
section. 

Theorem 3.1 Let T, xq, U be as in the further assumptions. Also let u G C‘^{U) 


satisfy 


\Lu\ < m U 

d>^u\r = 0 {\p\ < 1 ). 


( 3 . 1 ) 


Then there exists a small open neighborhood of U' <Z U of Xq so that n = 0 in U'. 

Next we will describe the Carleman estimate which we will use to prove Theorem 
o Let r]^ be the unit co-normal vector of T. By an appropriate translation and 
rotation, we can assume that xq = 0 and 77 ° = ( 1 , 0 , • • • , 0 )"'". We now make a change 
of coordinates near 0 by using the “Holmgren transform”, i.e., 

Xj ^ Xj, 2 < j < n, t = xi + k\x'\'^ for an appropriate constant k, 

where x' = {x2,--- ,Xn). In the new coordinates x = {x^,x') = {xi,X2,--- ,Xn) = 
it,x'), fl 3 .ip becomes 


\Lu\ < cJ 2 \,y\<i\Dxu\ in U 

supp(n) (Z {x : t > k|x'P} for some constants c, k > 0 , 


( 3 . 2 ) 


where u{x) = u{x{x)), 27 is a small neighborhood of 0 and the partial differential 
operator L is defined by 



l<p, q<n 1</9<W 
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with 


Cl'pix) = J \x{x)) C^Jf^{x{x))djXpdeXq, 

l<j, i<n 

the Jacobian J{x(x)) of the change of variables and partial derivative dp with respect 
to Xp. It can be easily checked that L satisfies the basic assumptions and further 
assumptions in U. 

Theorem 3.2 There exist positive constants To, ko, r, and c such that for 0 < T < 
Tq, we have for k >ko that 

^r-^(fcT)3-2hl /" < c /" \\Lvf dt (3.4) 

\u\<l -^0 -^0 

holds for all v{t,x') G with supp(n) C {{t,x') : |x'| < r, 0 < f < T/2}, 

where || • ||^ = (•, •) is the norm. 

Theorem 13.21 will be proved in the remaining sections. Once having this Carleman 
estimate, the proof of Theorem 13.11 is rather standard (see for instance [S]). We give 
the proof here for the sake of completeness. 

Proof of Theorem \3.1\ Let 6{t) G C°° defined in f > 0 with 6'(t) = 0 for t > T and 
9{f) = 1 for t < 2T/3. Assume that m is a solution of fl3.2p . When T is sufficiently 
small, we can apply fl3.4p to 9u and get that 

^l^l<,T-nA;T)3-2lH e^^^-'^^"\\d'^{9u)fdt < cj^ e^^^-^'>"\\L{9u)fdt 

< c (3.5) 

+c [L,9]u\\'^dt, 

where [L, 9] denotes the commutator of L and 9. By using the equation in (13. 2 p and 

taking k large, T small, if necessary, we can absorb the hrst term on the right side 

of (13.5p . Hence we have 

/•t/2 _ i>T 

Y / e'^^^-^^'\\d^ufdt<c'Y / e’^^^-^'^'Wd'^ufdt (3.6) 

for some P > 0. Since is decreasing in 0 < f < T, we have from fl3.6p that 

/•t/ 2 pT 

Y / / wm^dt, 

which implies that M = 0int<T/2. □ 
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4 Associated constant coefficients operators 


In this section we provide some factorization of L and consider an associated constant 
coefficient operator. To avoid any fnrther heavy notations, we snppress nsing in 
this sections bnt also in the rest of sections except the last section. That is we abnse 
the notations L etc. to denote L etc. 

From the basic assumptions, we have the factorization fl2.7p with simplified ex¬ 
pression on the dependency of rj, That is by taking r; = (1, 0, • • • , 0), ^ = (^i, ^') = 
(A,^') with = (^ 2 ; • • • ;Cn) ^ wc Can simplify the notations, for example M 

and Ph can be written as M = M{x,^) and Ph = Ph{x,^'), respectively. Since it is 
enongh to obtain a Carleman estimate for , we simply denote it by L. 

On this occasion, we introdnce the notations D = —i[d/{dxi), ■ ■ ■ ,d/{dxn)) and 
D' = —i[d/{dx 2 ),--- ,d/{dxn))- Fnrther we note here that we sometimes write t 
instead of Xi, that is t = Xi. 

Now we conjugate C := —L by the Carleman weight function Wk{t) := exp{k{t — 
r)^/2) with a large parameter k > 0. That is to consider 

Ck ■■= Wk{t)Cwk{t)~^. (4.1) 

The principal symbol of Ck under the scaling 

(A,^',fc) ^ (pA,p^',pfc) for p > 0 (4.2) 


is given by 

Hk = Hkix, e, k) = {x- - i{Z + kit - T))]}{A - [y + liZ - kit - T))]}. (4.3) 

We first fix xq € t/ and try to derive a Carleman estimate for C obtained from C 
by freezing its coefficients at xq. Note that C is given by 

£(A, D') = {A - iYiD'Y - iZiD'))}{Dt - (y(A) + y(A))}, (4.4) 

where y(f) = y(xo,C'), ^(^') = ^(a^o,f)- By dehning ij(0 = Hixo,^), we have 

Ak = AkiDt, D') = Wkit)HiDt, D’)wkit)~'^ = Pk,bPk,g, (4.5) 

where the good part Pk,g and bad part Pk,b are given by 

Pk,g = A - GkiD’), Pk,b = A - GliD') (4.6) 

with 

GtK') = Y((') + - k(t - T)), G’M') = F(ey - + k(t - T)). (4.7) 

To have an a priori estimate for ^ and Pk^b, we construct right parametrices 
Sk,g and Sk,b for their adjoint operators P^^ and P^ respectively. Since P^^ is an 
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elliptic operator of order one with large parameter fc, Sk^g can be easily constructed 
using the theory of pseudo-differential operator with large parameter to have 


Pk,g^k,g — I + Rk,g (4.8) 

for some pseudo-differential operator R^ g of large negative order with large parameter 
k which satisfies the estimate 

\\{I + Rk,g)-^\\<1 {kT>Mg) (4.9) 

for some Mg > 0, where the norm || • || is the operator norm on Further, for 

each s G M, there exists a constant Cg(s) > 0 such that Sk,g satisfies the estimate 

||‘5A:,c,'?n||H'*(R") < Cg{s){kT) {w G if^(M”)). (4.10) 

As for Sk,b we construct a right pararametrix S'f. ;= ^ for 

n = Kb = ^Pkfi (4-11) 


in the next section. 


5 A right parametrix of VI 

We look for the Schwartz kernel *S'^(xi, yi; ^') of a right parametrix iS(, of the ordinary 
differential operator in xi with parameters k near its characteristics in the 
form 

S't(x„y,X) = (5.1) 

where 

4>(a;i,2/i,f) = (xi-2/i)^i + $i(a:i,|/i,f) 

=7G;7(«i-y(n->(z(5')-Mr-9i))) 

and d^i = (27r)“^d.^i, rj G C^(M), r]{s) = 1 for |s| < 1/4 and r]{s) = 0 for |s| > 1/2 
with 0 < ?7 < 1. Further 7 > 0 is defined as follows. The principal symbol of o'{P^) 
of P/ is given by 

a(Fi) =?,-r(r)-i(7e')-Mr-i,))) 

= - {B((') - tk(T - ^' ’ 

where B = B{xo,^') is the simplihed notation of B{xo,ri,^) given by fl2.8p for rj = 

(l,0,---,0),e=(ei,n = (A,f). 

Due to the symmetry of Br{x, ^') and Bj{x, ^'), observe that for any v G Ker cr(P/) 
with |n| = 1, we have 

iiiSjia-HT-Xi))v,v) = 0 (5.4) 
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and 


(te-5ii(nK^) = 0. (5.5) 

Further since Bi{x,^') = Z{x,^') is positive, we have for xi G [0,T/2], 

X-\kT) < le'l < X{kT) (5.6) 

for some constant A > 0 from fl5.4p and 

161 < ”^161 with m = max{|in 6 |=i(.BHn,n)|, |sup|^|=i(^i?^^) ^^)l} (5-7) 

from 05.51) . Combining 05.7p with 05.6p . we have 

161 6 l{kT) with 7 = Am. (5.8) 

Note that m can be taken uniformly for xq G U. 

Based on these the characteristics of is simply given by 

det{6 - (5(^0 - tk{T - xi))} = 0 (5.9) 


for large k > 0 and it is enough to consider the operator with prinicipal symbol 
6 - (5(6) - ik{T - a:i)) and (T - |/i )"^(6 - (5(6) - *^(5 - yi)) instead of VI 
and respectively. In the rest of this paper, we denote them by and 

^(l/i) 6 )) respectively. 

Since we want to have 


^T—xi 

we can take 


Now let 


= [iirKM. + At - + ik]M-'(!/i,a 


T—x\ 


T—x\' 




(5.10) 

(5.11) 


S ' k , i { x , y ) = /r.-i ^a(| 5 ) 6 (^' 62/116)^^6)^ .g^2) 

= / r.-i s'l^ixi, 2 / 1 ; 6)C 


and 




(5.13) 


where d^' = (27r) ^6^, d^ = (27r) and G C^(M"' ^) such that 0\{^') = 

1 ((2A)-1 < 161 < 2A), 066) = 0 (6 ~ 0) and 60 = 1-6 with 6(6 = e^{i'/{kT))y{i,/{2^kT)). 
Then, by decomposing S'^{x,y) into S'i.{x,y) = 6 1 (^ 1 2/) + 6o(^i2/)i we will prove 
that 6 (^ 12 /) is the Schwartz kernel of a right parametrix 6 of VI in the next section. 

Before closing this section, we remark here that we have 


|det6(e)V(!/i.?'))l >r-''(|{| +1;T)", (5.14) 

where the notation denotes ”>” modulo multiplication by a general positive 
constant. 


14 














6 Error estimate of parametrix S[ 


k 


In order to see that S'/, is a parametrix of we need to estimate the error 

n‘4 - 1 - (6-1) 

Since / = 6{x — y) with d^ = {2Tr)~^d^ and taking 

into account fld.lOj) . the sum of two errors coming from the operators with Schwartz’s 
kernels Sl^{x,y), S^^ix^y) is 


+ 4^-1 {J “ v) 

in terms of Schwartz kernel. 

Next we will show that by taking k and T to satisfy k > T~^, we have 

= - !/i) + 0(fc-‘/») 


( 6 . 2 ) 


(6.3) 


over a neighborhood of Xi = yi uniformly respect to any hxed G yi G 

[0,T/2], where 0{k~^/^) is the kernel of a bounded operator on L^((0,T/2)) norm 
with operator estimated by 0((fc)“^/®) uniformly with respect to other variables such 
as 1 / 1 , For this it is enough to show 




(6.4) 


over a neighborhood of Xi = //i, where 4 * 2 ( 2^15 l/i, Ci) = [xi — yi)ii — ^ 2 (T-yl) '^1 

43(xi,!/,,e) = f^BK'). Write 

*2(a:i,9i,?i) = (ii -9 i)5i - = (^i - with Wi = (1 - ^^E|i))?i. 

(6.5) 

Since the Jacobian of the transformation uji hg- is (1 — ^44)) ) have 

= «a(|.) /e‘<"--«’“‘>>)((27fcr)-‘(l - 

( 6 . 6 ) 

This is the Schwartz kernel of the pseudo-differential operator over (0,T/2) with 
double symbol 

z := »A(|r)ii((27*:T)->(l - 5^)-‘iei)(l - (6.7) 

and its simplihed symbol Zl has an asymptotic expansion 

ZL{x,bJi,C) ~ 

E.(a!)-’a2;r'?;{7((27fcr)-Hi - 53e|^)-V)(i - 5£^)-ie»3ta-...5')} . 

( 6 . 8 ) 
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Here note that gives an effect of amplifying by {kT) ^ and that of Dy^ is ampli¬ 
fying by max{fc^/^, T“^}. These total effect can be transformed into amplifying by 
0 ((fci/ 2 T)-i) if we take k and T to satisfy k >T 
Thus, we have 

V:s', = l + R,, (6.9) 

where Rb is a bounded linear operator on T^((0, T/2) x with its operator norm 

estimated by 0{k~^^^) if we take k and T to satisfy k > T~^. Therefore we can modify 
iS(, to obtain a fundamental solution of V^. which can have the same estimate as that 
of 5' if we can show that is a bounded linear operator on L^((0,T/2) x M"' ^). 
Also, this relation on the estimates of iS(, and its fundamental solution is true for 
their derivatives. 


7 Estimate of (xi, 


Let xo G 17 be such that B(^') = B(xo,(') is diagonalizable for any G M” ^ \ 0. 
Put Bb = B + ik{yi — T). Then 

M(!/,= (T - !„)(6 - B,)-\ = -|HirT(6 - B). 

By the further assumptions and spectral decomposition of H, we have 




(7.1) 


where Xh = \h + ik{yi — T), Ph = Ph{xo,^') is the projection defined before by (12.111) 
and (p(A) is defined by 


(p(A) = exp[i(xi - yifX/ (2(T - |/i))]. (7.2) 


Let \j = A(xo, CO = Vj+ijUj with Vj = CO; hj = hji^o, CO G M for 1 < j < s 
be the eigenvalues of H(C 0 = B{xo, CO with multiplicity rrij satisfying mi + - ■ ■ + ms = 
N. For the estimate of fl7.ip . we need to estimate 






^ ICl “ — Rhj 


A)(T-|/i))}-M~Ci. 


By taking 7 large to satisfy |z/0 < l/^i'jkT) (1 < j < N) so that each of |(Ci + 
Uj)/{2'jkT)\ < 1/2 (1 < j < iV) implies \^i/{2'jkT)\ < 1, it is enough to estimate 


nj(xi,|/i;xo,C0 := 


fc(T-2/i))}-MCi 
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for each 1 < j < iV. For the estimate of each observe that 


I I\^i\<2'ykT 
<(27r)-V. 

= 7r~^ log 2 . 




'■ 2'ykT ^ 
lkT<\ii\<2"ikT 


So we only need to estimate 


'^2'ykT / nz 

^i{{xi-yi )-——^ 

J —2'ykT 

To begin estimating this, set 


gd(^i yi) 2 ?Td;V)Fie - A;(T - 


(7.3) 


I = %-yl} 1^1 e 2 (rT,Y) (7.4) 


(^1—yi) 


/ — 2'ykT 


and y = yj{yi, x = Xi - yi, m = T - yi, ( = ^i, dC = d^i, then 


I = /! yiK^-x^m V/2 . _ ^(^ _ km)}-^dC 

_ p—x^mr"^P l/ 2 r2'ykT i(x—x‘^m~^ /2)C C~ki{lJ:—km) 

J-2'rkT^ C^ + {^i-kmf^^^ 


Now, we consider two cases. For the case y — km 7 ^ 0. Observe that 
\y - km\ = tan-H^^)]?: 

So, we only need to estimate 

J ._ 1 r2'fkT 


2'ykT — 


(7.5) 


(7.6) 


J. J. PV' < i(a.-»>m-V2)C Ji-.XrO 

2i ' 2i Jo fcm)^ ’ ^ 2i J — 

_ r2'ykT Q-s\n[{x—x‘^m~^/2)Q j^ 

Jo C,'‘^ + {yt—kmy^ ’ 

_ r2'rkT C-sin[bC] jy 

- Jo ( 2+^2 

_ j2ybkT (sinC 


__ ( __ pi{x-x'^m V2)C y/- 

2'ykT C2+(^_fcm)2 


^0 


(^+b'^a^ 


dC, 


(7.7) 

where a = y — km, b = x — x^m ^/ 2 . 

If 6 = 0, then we just have J j {2i) = 0. Hence we assume & 7 ^ 0. Let / G N U {0} 
such that 2'ybKT — 2nl < 27r. For any natural number j [j < /), we have 

I (’27rj+27r (.sint^ j/-\ 

\ J2nj f2+b2a2“SI 

r27rj+7r (^.sinC i r27rj+27r (.ginC J AI 

J27Tj C2+fe2a2“S J2^j+,r (2 _pi,2^2 \ 

r27rj+7r TsinC 7^ I f27rj+7r (i^+^).sin(C+7r) | 

J2'7Tj C‘2-\-b2a2 ^ 72771 ('C-l-7r')2-l-(>2a2 > I 

r2'Kj 

J2'Kj 


27rj+7r C-sinC yy _ 

^ J27r7 


C^+fe2a2 

r27rj+7rr_^ 

72777 1 


277J ((^+77)24.f)2 a2 

277J+77 (C+77).sinC 


27''J ((^+77)2+b^a2 

—'.sinCdC 


(7.8) 


)2'nj l.C 2 +fe 2 (j 2 (^+,r)2+fe2 

r277j + 77 C^77 + C77^-fe^a 77 • 

7277; (C2+fe2a2)[((^+77)2+b2a2] 

< TT^dC. 

- J 27TJ lOd-TTl-^ J 


(C+77)2 
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Now, we estimate We have that 


|Z| — |V^-1 r27rj+27r ^-sinC , r^ir (sinri , r 
I 2il “ IW=1 J27rj (2+62^2 Jo (2+62„2«t, -t- Ja^z 

<rcllF<iC + 4>r 

< 2 + 47r. 


2’ybkT ^-sin^ 


C2+f,2c 


(7.9) 


On the other hand, for the case /x — km = 0, we have 


j| ^ V/2 


_ |2 g-a:^m“^/4/2 sm[{x-x^m ^)C] 


r 

Jo 


C 


idCI- 


The same argnment can give |/| bonnded by an explicit positive constant. 


{^x\—y\) (c/\ 

Since (1 + k[xi — ^ e~ for m > 0, we have hnally obtained 

the estimate for S'^{xi,yi]i') = S[{xi,yi]XQ,i') 


(7.10) 


\s',{xr,y,-,0 \<{^ + K^i-yiYr"T 


(7.11) 


uniformly for any Xi, yi G [0,T/2], G \ 0 and of those Xq E U with the 

property stated in the first line of this section. Here T in the above estimate comes 
from the factor T — yi in fl7.ip . 


8 Estimate of parametrix 

Recall that the Schwartz kernel of iS(, was given by 


S'ki.x,y) 


=/m- 1 ^A(|^)e'("' 

= e^^^'-y'>^'Sk(xi,yi;C)dC- 


) ,yi .«') M-i (yi , e') ) d^' 


In this section, we will estimate d^S'^^v for |q;| < 1 and g G O“((0,T/2) x R"" ^). To 
begin with, 

\\Skg\\L^iTl^-\dx') = (/|*S^c/(x)|W)^/2 

= (/ I II Sk{x,y)g{y)dy'dyi\^dx'y/^ (8.1) 

<IiI\I S'kix, y)g{y)dyfdx'Y/^dyi, 

where we used the Minkowski integral inequality in the last inequality. 

Now we estimate the inner integral 


I S'k{x,y)g{y)dy' 

= II yi, i')g{y)dy'di' 

= IjIY"'^'S'^{ xi,yi,i')g{yi,i')di' 

= S'kg{yi,x), 
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where ” • ” denotes the inverse Fonrier transform of 
and fl7.1ip . we can estimate 


By the Plancherel theorem 


So, we have 


(/1 1 y)9{y)dyVdx'y/^ 

< T[1 + k{xi - yif]-^ • (/ \g{y)Wy/^. 


ll‘^fc5'llL2((o,T/2)xR"-i) ~ fo ll‘^A:5'(^)llL2(Rn-i 

< / [ / ^[1 + ^(^1 - yi?]~^ ■ (I \9{y)?dy'f/‘^ dyi] ^dxi 

< T\J[1 + kx‘l]-^dxif\\g\\l2^j^n) 

2 

L2((o,r/2)xR"-i)’ 


(8.3) 


(8.4) 


< V 


where we used the Young inequality in the second inequality. 

Let ^ and 5^ q be the operators with Schwartz’s kernels y) and g(x, y), 

respectively. Then, the same argument gives that for |q!| < 1 

ll^“‘5fc,lfi'll^h(0,T/2)xR—1) < :^(^^)'“'llfi'l|L2{(o,T/2)xR"-i)- (8.5) 

For the estimate of qU for |q;| < 1, we surely have from fl5.14p that 

ll^“‘5fc,o5'l|L2{(o,T/2)xR"-i) < :^(^^)'“'||5 'I|l2((0,T/2)xR"-1)- (8.6) 

Thus for I a I < 1, we have 

ll^“‘5fc5'IU2((0,T/2)xR"-i) < :^(^^)'“'||5'IU2((0,T/2 )xR"-1)- (8.7) 

Combining fl8.7p with what was mentioned in the last paragraph of Section 6, we 
obtain for v e C'^((0,T/2) x M"'“^) that 




*"-i) ^ ^ll”Pfc^'l|L2((0,r/2)> 




~ l|-Pfc,fe'^^IU2((o,T/2)xR’>-i)- 


( 8 . 8 ) 


Since can have a similar estimate as fl8.7p and its adjoint is (—l)l"l9"S'^, we can 
see likewise fl8.8p 


X]|/3|<1^ 2(/cT) 2 ||9fn||L2((o,T/2)xR'i-i) ^ ||-Pfc,6^^||L2((o,T/2)> 


(8.9) 


9 Estimate of A 


k 


First of all since Pk^g is an elliptic pseudo-differential operator of order 1 with large 
parameter k, we can easily have 








for any v G C^((0, T/2) x M"" by the standard argnment (see for instance Chapter 
2 of [^). 

Combining fl8.9p and flh.ip . we have for any v G C^((0, '^) x that 

< \\A^vf ( 9 . 2 ) 

and hence we have the following Carleman estimates for constant coefficients partial 
differential operator C. 

Theorem 9.1 There exist positive constants Tq, ko and c such that for 0 < T <Tq, 
we have for k > ko that 

{kT‘^)~^ + k e'^^^^~'^^^\\dv\\‘^dxi + k^T'^ e'^^^^~'^^^\\v\\‘^dxi 

< c e^^^'^~'^'^^\\Cv\\‘^dxi 

(9.3) 

for any V G C'^((0,T/2) x M"' where || • ||^ = (•, •) is the L‘^{MT norm. 


10 Carleman estimates by partition of unity 

We will nse a partition of nnity and Theorem 19.11 to prove Theorem 13.21 To begin 
with we first note that we are back to nse the notations nsed to describe fl3.3p . let 
■do G C^(M) and 0 < "d < 1 snch that 


do(t) 


1 (|t| < 1) 

0 (|t|>3/2). 


Also let {}{x) = 'do(a;i) • • •'do(a;„). Then, we have 


■d^x) 


1 (a; G Qi(0)) 

0 (5iGM'^\g3/2(0)), 


( 10 . 1 ) 


( 10 . 2 ) 


where Qr = gr(0) = {x G M” : \xj\ < r, j = 1,2, ■ ,n}. Farther for /r > 1 and 

(7 G Z”, dehne 

Xg = g/ti 

and set 

^gA^) = - g)- 

Then we have 

snppdg,^ C Q3/2g{Xg) C Q2/g{Xg) (10.3) 

and 

1-^ '^g,r-\ — (XQ3/2^(ig) ~ XQi/^{xg)) {k = 0,1, 2), (10.4) 

where Ci > 1 depends only on n. 
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For G Z"', let 


Ag = {g' e Z*" : supp'd n supp'dg^g 0}, 
then the number of Ag will only depend on n but not on p. Thus, we can define 

:= T^gez^^gA^) > 1 e M’"). (10.5) 

From fllO.dp . we have 

< C2fi\ (10.6) 

where denotes all the k th order derivatives and C 2 > 1 depends only on n. 
Dehne 

VgAA = ^gAAAgA) (x e M”). 

Then, we have 

f Egez^ VgAA = 1 (5 e M^) 

< suppr/g,^ C Q3/2t,(Xg) C Q2/i,(Xg) (10.7) 

[ \DAg,g\ < C3/XQ3/2Mfe) = 0, 1 , 2 ), 
where C 3 > 1 depends only on n. 

Now define Lg{x,d)u = 'E.i<j,£<n^i<p<N Then from Theorem 

19.11 we have the following proposition. 

Proposition 10.1 There exist positive constants Tq, ko, r and c such that for 0 < 
T <Tq, we have for k > ko that 

{kT^)~^ e’^^^^~'^'>^\\d'^v\\‘^dxi + k e^^^^~'^'>^\\dv\\‘^dxi + k^T"^ e^^^^~'^'>^\\v\\'^dxi 

< c e^^^^~'^^^\\Lgv\\‘^dxi 

( 10 . 8 ) 

for all v{xi,x') G (^^(M") with suppn <ZU ■.= {(xi,x') : Xi G (0,T/2), \x'\ < r} C 

U. 

Now let supp(n) C W/2 and apply fllO.Sp to vpg^g. Then we have 

(/cT2)"^/(f + k e^^^^-'^'^^\\D(ygg^g)\\'^dxi 

< cj^ e^^^^-^^^\\Lg{vgg^g)fdxi. 

Note here we have 

\Lg{x,d){ygg^g)\< \L{x,d){ygg^g)\ 

+ \L{x,d){vgg,g) - Lg{i,d){vgg^g)\ 

< ggjL{T,d){v)\ + pADMxq2{x,) 

+t\Dv\ Xq2{x,) + A\v\ Xq2{xs)- 


(10.9) 


( 10 . 10 ) 
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By summing up all fllO.lOp with respect to (7 G Z” and using from fll0.7p . fll0.9p . fll0.10p . 

we have 

[(/cT^)“^ — + [k — C2/i^] e^^^^~'^^^\\Dv\\^dxi 

+ [k^T^ - c^kfil-cfi^] 

( 10 . 11 ) 

Taking C 2 p^ = k/2, AT = c^jcx and k large, we can derive that 

{kT‘^)~^ \\D‘^v\\‘^dxi + k WDvW^dxi 

+/c^T^/g e'^^^^~'^^^\\v\\‘^dxi ( 10 . 12 ) 
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